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Critical behavior of repulsive dimers on square lattices at 2/3 monolayer coverage 
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Monte Carlo simulations and finite-size scaling theory have been used to study the critical behavior 
of repulsive dimers on square lattices at 2/3 monolayer coverage. A "zig-zag" (ZZ) ordered phase, 
characterized by domains of parallel ZZ strips oriented at ±45° from the lattice symmetry axes, 
was found. This ordered phase is separated from the disordered state by a order- disorder phase 
transition occurring at a finite critical temperature. Based on the strong axial anisotropy of the ZZ 
phase, an orientational order parameter has been introduced. All the critical quantities have been 
obtained. The set of critical exponents suggests that the system belongs to a new universality class. 
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I. INTRODUCTION 

The study of critical phenomena and phase tran- 
sitions is a major and long standing topic in sta- 
tistical physics. Particularly, the two- 
dimensional lattice-gas model 8 - with repulsive interac- 
tions between the adparticles has received considerable 
theoretical and experimental interest because it pro- 
vides the theoretical framework to study structural order- 
disorder transitions occurring in many adsorbed mono- 
layer £1^^11,1^3,14,15,16,17,18,19,20,21,22 Most studieg 

have been devoted to adsorption of particles with single 
occupancy. The problem becomes considerably difficult 
when particles occupy two adjacent lattice sites (dimers). 
Consequently, there have been a few studies devoted to 
order-disorder transitions associated to dimer adsorption 
with repulsive lateral interactions. Among them, the 
structural ordering of interacting dimers has been ana- 
lyzed by A. J. Phares et al^ 3 - The authors calculated the 
entropy of dimer on semi-infinite M x N square lattice 
(N — > oo) by means of transfer matrix techniques. They 
concluded that there are a finite number of ordered struc- 
tures. As it arose from simulation analysis r 4 - only two of 
the predicted structures survive at thermodynamic limit. 
In fact, in Ref. [24J, the analysis of the phase diagram for 
repulsive nearest-neighbor interactions on a square lattice 
confirmed the presence of two well-defined structures: a 
c(4 x 2) ordered phase at 9 = 1/2 and a "zig-zag" (ZZ) 
order at 9 — 2/3, being 9 the surface coverage. 

The thermodynamic implication of such a structural 
ordering was demonstrated through the analysis of ad- 
sorption isotherms^ the collective diffusion coefficient^ 
and the configurational entropy^ of dimers with nearest- 
neighbor repulsion. Later, Monte Carlo (MC) simula- 
tions and finite-size scaling (FSS) techniques have been 
used to study the critical behavior of repulsive linear 
fc-mers in the low-coverage ordered structure (at 9 = 



l/2) i 28 ' 29 A (2k x 2) ordered phase, characterized by al- 
ternating lines, each one being a sequence of adsorbed 
fc-mers separated by k adjacent empty sites, was found. 
The critical temperature and critical exponents were cal- 
culated. The results revealed that the system does not 
belong to the universality class of the two-dimensional 
Ising model. The study was extended to triangular lat- 
tices^ In this case, the exponents obtained for k > 1 
and 9 = k/(2k + 1) are very close to those characterizing 
the critical behavior of fc-mers (k > 1) on square lattices 
at 9 = 1/2. 

Recently, by using MC simulations and finite- 
scaling techniques, Rzysko and Borowko have stud- 
ied a wide variety of systems in presence of multisitc- 
occupancyi 31 ' 32 i 33 ' 34 i 35 Among them, attracting dimers 
in the presence of energetic heterogeneity; 3 -! heteronu- 
clear dimers consisting of different segments A and B 
adsorbed on square lattices ; 32 ! 33 ' 34 and trimers with dif- 
ferent structures adsorbed on square lattices^ In these 
leading papers, a rich variety of phase transitions was 
reported along with a detailed discussion about critical 
exponents and universality class. 

Summarizing, although there have been various studies 
for monolayers at half coverage, to the author's knowl- 
edge, there are no conclusive studies on the characteris- 
tics of the transition phase of repulsive dimers on a square 
lattice at 2/3 coverage. In the present contribution we 
attempt to remedy this situation. For this purpose, ex- 
tensive MC simulations in the canonical ensemble com- 
plemented by analysis using FSS techniques have been 
applied. The FSS study has been divided in two parts. 
Namely, 1) a conventional FSS in terms of the normal- 
ized scaling variable e = T/T c — where T c is 
the critical temperature; and 2) an extended FSS, 38 i 39 
where a = 1 — T c /T, instead of e, is used. Our results led 
the determination of the critical temperature separating 
the transition between a disordered state and the ZZ or- 
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FIG. 1: Snapshot of the ordered phase for dimers at 9 = 2/3. 



dered phase occurring at 2/3 coverage and the critical 
exponents characterizing the phase transition. 

The outline of the paper is as follows: In Sec. II we 
describe the dimer lattice-gas model. The order param- 
eter and the simulation scheme are introduced in Sees. 
Ill and IV, respectively. Finally, the results and general 
conclusions are presented in Sec. V. 



II. THE MODEL 

In this section, the lattice-gas model for dimer adsorp- 
tion is described. The surface is represented as a simple 
square lattice in two-dimensions consisting of M — Lx L 
adsorptive sites, where L is the size of the system along 
each axis. The homonuclear dimer is modelled as 2 
monomers at a fixed separation, which equals the lat- 
tice constant a. In the adsorption process, it is assumed 
that each monomer occupies a single adsorption site and 
the admolecules adsorb or desorb as one unit, neglecting 
any possible dissociation. The high-frequency stretching 
motion along the molecular bond has not been considered 
here. 

In order to describe the system of N dimers adsorbed 
on M sites at a given temperature T, let us introduce 
the occupation variable Cj which can take the following 
values: Cj = if the corresponding site is empty and 
Cj = 1 if the site is occupied. Under this consideration, 
the Hamiltonian of the system is given by, 



H 



Nw - 



(i) 



where w is the nearest-neighbor (NN) interaction con- 
stant which is assumed to be repulsive (positive), (i, j) 
represents pairs of NN sites and e is the energy of ad- 
sorption of one given surface site. The term Nw is sub- 
tracted in Eq. (|TJ) since the summation over all the pairs 
of NN sites overestimates the total energy by including 
N bonds belonging to the N adsorbed dimers. 



III. ORDER PARAMETER 

Given the inherent anisotropy of the adparticles, it is 
convenient to define a related order parameter. In this 
section, we will briefly refer to a recently reported or- 
der parameter <5 j22 which measures the orientation of the 
admolecules in the ordered structure. 

Fig. Q] shows one of the possible configurations of the 
ordered ZZ structure appearing for dimers at 2/3 mono- 
layer. Though the degeneracy of this phase is high, the 
entropy per lattice site tends to zero in the thermody- 
namic limit^I. The figure suggests a simple way to build 
an order parameter. In fact, any realization of the ZZ 
structure implies the orientation of the particles along 
one of the lattice axis4£ Then, all the available config- 
urations can be grouped in two sets, according to this 
orientation. Taking advantage of this property, we define 
the order parameter as: 



S = 



N 



(2) 



where N v (Nh) represents the number of dimers aligned 
along the vertical (horizontal) axis and N = N v + Nh- 

When the system is disordered (T > T c ), the two ori- 
entations (vertical or horizontal) are equivalent and S 
is zero. As the temperature is decreased below T c , the 
dimers align along one direction and 5 is different from 
zero. Thus, 5 appears as a proper order parameter to 
elucidate the phase transition. 



IV. MONTE CARLO METHOD 

The lattices were generated fulfilling the following con- 
ditions: 

1) The sites were arranged in a square lattice of side L 
(M = L x L) , with conventional periodic boundary 
conditions. 

2) Because the surface was assumed to be homoge- 
neous, the interaction energy between the adsorbed 
dimer and the atoms of the substrate e Q was ne- 
glected for sake of simplicity. 

3) In order to maintain the lattice at 2/3 coverage, 
8 = 2N/M = 2/3, the number of dimers on the 
lattice was fixed as N = M/3. 

4) Appropriate values of L (= 60, 72, 84, 96, 108) were 
used in such a way that the ZZ adlayer structure is 
not altered by boundary conditions. 

In order to study the critical behavior of the system, 
we have used an exchange MC methodj 41 ' 42 As in Ref. 
EL we build a compound system that consists of m non- 
interacting replicas of the system concerned. The z-th 
replica is associated with a heat bath at temperature 
[or [3i = 1/(/cbTj), fcs being the Boltzmann constant]. 
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To determine the set of temperatures, {T}, we set the 
highest temperature, T±, in the high-temperature phase 
where relaxation (correlation) time is expected to be very 
short and there exists only one minimum in the free en- 
ergy space. On the other hand, the lowest temperature, 
T m , is set in the low-temperature phase whose properties 
we are interested in. Finally, the difference between two 
consecutive temperatures, T and T + i with Tj > T + i, is 
set as AT = (T\ — T m ) j(m — 1) (equally spaced temper- 
atures) . 

Under these conditions, the algorithm to carry out the 
simulation process is built on the basis of two major sub- 
routines: replica-update and exchange. 
Replica-update: Interchange vacancy-particle and diffu- 
sional relaxation. The procedure is as follows: (a) One 
out of the m replicas is randomly selected (for example 
the i-th replica), (b) A dimer and a pair of nearest- 
neighbor empty sites, both belonging to the replica cho- 
sen in (a), are randomly selected and their coordinates 
are established. Then, an attempt is made to inter- 
change its occupancy state with probability given by the 
Metropolis rule, — : 



P = mm{l,exp(-ftA/T)} 



(3) 



where AH is the difference between the Hamiltonians 
of the final and initial states, (c) A dimer is randomly 
selected. Then, a displacement is attempted (following 
the Metropolis scheme), by either jumps along the dimer 
axis or reptation through a 90° rotation of the dimer axis, 
where one of the dimer centers remains in its position 
(interested readers are referred to Fig. 1 in Ref. [26| for 
a more complete description of the reptation mechanism) . 
This procedure (diffusional relaxation) must be allowed 
in order to reach equilibrium in a reasonable time. 
Exchange: Exchange of two configurations Xi and Xy, 
corresponding to the i-th and i'-th replicas, respectively, 
is tried and accepted with probability W (Xi, ft|-XV , ft'). 
In general, the probability of exchanging configurations 
of the i-th. and i'-th replicas is given by, — 



W (Xi, (3i\Xi> , Pi>) 



1 

exp(- 



for 
-A) for 



A < 
A > 



(4) 



where A = (ft - ft,) [#p£» - ff(Xj)]. As in Ref. 
we restrict the replica-exchange to the case i' = i + 1. 
The complete simulation procedure is the following: 

1) Initialization. 

2) Replica-update. 

3) Exchange. 

4) Repeat from step 2) m x M times. This is the 
elementary step in the simulation process or Monte 
Carlo step (MCS). 

The initialization of the compound system of m repli- 
cas, step 1), is as follows. By starting with a random 



TABLE I: Parameters used in the MC runs. For all lattice 
sizes we have chosen T\ = 0.25 and T m — 0.15 (the tempera- 
tures are in units of w/ks)- 



L 


m 


ni 


U2 


UMCS 


60 


101 


10 5 


10 b 


3 x 10 b 


72 


121 


10 5 


10 5 


3 x 10 s 


84 


141 


2 x 10 5 


2 x 10 s 


5 x 10 5 


96 


161 


4 x 10 5 


4 x 10 s 


8 x 10 5 


108 


141 


5 x 10 5 


5 x 10 s 


10 6 



initial condition, the configuration of the replica 1 is ob- 
tained after m MCS' at T x (MCS' consists of M real- 
izations of the replica-update subroutine). Second, for 
i = {2, m}, the configuration of the i-th replica is ob- 
tained after ri\ MCS' at T, taking as initial condition 
the configuration of the replica to T%—\. This method re- 
sults more efficient than a random initialization of each 
replica. 

The procedure l)-4) is repeated for all lattice sizes. For 
each lattice, the equilibrium state can be well reproduced 
after discarding the first rii MCS. Then, averages are 
taken over umcs successive MCS. As it was mentioned 
above, a set of equally spaced temperatures is chosen in 
order to accurately calculate the physical observables in 
the close vicinity of T c . 

The thermal average (...) of a physical quantity A is 
obtained through simple averages: 



(A) 



1 



TlMCS 



riMCS 

*=i 



(5) 



In the last equation, Xi stands for the state of the z-th 
replica (at temperature T). Thus, the specific heat C (in 
k B units) is sampled from energy fluctuations: 



1 



C 



k B (Lk B Tf 



l(H 2 ) (H) 2 ]. 



(6) 



The quantities related with the order parameter, such as 
the susceptibility x-> an d the reduced fourth-order cumu- 
lant U introduced by Binder, — are calculated as: 



k B T 



K<5 2 > - 



and 



17 = 1- 



(S 4 



3 (£ 2 



(7) 



(8) 



Finally, in order to discuss the nature of the phase 
transition, the fourth-order energy cumulant, Ue, was 
obtained as: 



U E = 1 



(H 4 



(9) 



4 
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FIG. 2: (Color online) Size dependence of the order parameter 
as a function of temperature. Inset: dependence of (5) on 
i" 1 for different regimes of T as indicated. The error in each 
measurement is smaller that the size of the symbols. 

V. RESULTS AND CONCLUSIONS 

The critical behavior of the present model has been 
investigated by means of the computational scheme de- 
scribed in the previous section and FSS techniques. The 
values of the parameters used in each MC run are shown 
in Table I. In addition, all the simulation calculations 
were obtained by averaging over 50 MC runs. 

Because the replica temperatures were chosen equally 
spaced, the acceptance probability of the replica- 
exchange decreases in the critical temperature region, 
reaching a minimum whose value is always greater than 
50%. The equilibration has been tested by studying how 
the results vary when the simulation times n-i and umcs 
are successively increased by factors of 2. We require 
that the last three results for all observables agree within 
error bars. This simple method is shown to be useful to 
test equilibration [see, for instance, Ref. HJ. All calcula- 
tions were carried out using the parallel cluster BACO of 
Universidad National de San Luis, Argentina. This fa- 
cility consists of 60 PCs each with a 3.0 GHz Pentium-4 
processor. 

The conventional FSS implies the following behavior 
of C, (S), x an d U at criticality, 

C = L a/v C{L 1/u e) (10) 
(5) = L-^ u 5(L^ u e) (11) 
x = L^x(L 1/u e) (12) 

U = l>(L 1/ly e) (13) 

for L — > oo, e — > such that L l l v e= finite, where 
(e = T/T c — 1). Here a, ft, 7 and v are the standard 



FIG. 3: (Color online) Size dependence of the susceptibility 
as a function of temperature. The error in each measurement 
is smaller that the size of the symbols. 




FIG. 4: (Color online) U versus fcsT/w, for different sizes. 
From their intersections one obtained feTc/w. In the inset, 
the data are plotted over a wider range of temperatures. 
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FIG. 5: (Color online) Temperature variation of Ue for vari- 
ous lattice sizes. 
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FIG. 6: (Color online) (a) Log-log plot of the size dependence 
of the maximum values of derivatives of various thermody- 
namic quantities used to determine v. (b) Log-log plot of 
the size dependence of the maximum value of the suscep- 
tibility, the point of inflection of the order parameter and 
the maximum value of the derivative of the order parame- 
ter used to determine 7 and j3, respectively, (c) K C (L) vs. 
L _1// " from several quantities as indicated. From extrapola- 
tion one obtains the estimation of the critical temperature. 
In all cases, dotted lines correspond to linear fits of the data 
and L = 60, 72, 84, 96, 108. 



critical exponents of the specific heat ( C ~ |e|~ a for 
e — ► 0, L — > 00), order parameter ((<$) ~ —e" for e — > CP, 
L — > 00), susceptibility (x ~ |e| 7 for e — > 0, L — > 00) 
and correlation length £ (£ ~ |e| ~ w for e — > 0,i — » 00), 
respectively. (7, <5, x and U are scaling functions for the 
respective quantities. In the case of extended FSS, 38 i 39 
e and L in eqs. Ijl0lll3[) are replaced by a = 1 — T c /T and 
LVT, respectively. 
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FIG. 7: (Color online) Conventional data collapsing for: (a) 
the curves in Fig. [21 (b) the curves in Fig. [3] and (c) the 
curves in Fig. [4] 



We start with the calculation of the order parameter 
(Fig. [2]), susceptibility (Fig. [3]) and cumulant (Fig. QJ 
plotted versus ksT/w for several lattice sizes. Due to 
computational limitations^ the curves in Fig. [5] do not 
clearly show the existence, at thermodynamic limit, of 
a finite temperature below which the order parameter is 
different from zero. In order to clarify this point, the 
inset in Fig. [2] shows the dependence of (6) on L^ 1 for 
constant T. As it can be observed, (5) tends to a finite 
value (zero) for T < T c (T > T c ). 

From the intersections of the curves in Fig. 3] one gets 
the estimation of the critical temperature. In this case, 
ksTc/u] — 0.182(1), which is in good agreement with the 
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FIG. 8: (Color online) Same as Fig.|6]for the extended scaling 
scheme. 



FIG. 9: (Color online) Same as Fig.[7]for the extended scaling 
scheme. 



value previously reported in the literature^. In Ref. I24I 
the critical temperature was obtained from the peaks of 
the curves of the specific heat versus temperature (at 
fixed coverage) and coverage (at fixed temperature). A 
more rigorous study was not possible due to the lack of 
an adequate order parameter. In the inset, the data are 
plotted over a wider range of temperatures, exhibiting the 
typical behavior of the cumulants in presence of a contin- 
uous phase transition. With respect to the value of the 
cumulant at the transition temperature, U* , this quan- 
tity was calculated by plotting U*(L) vs L" 1 ^—, where 
the value of U*(L) was obtained by fixing Ki= w/hgT) 
at our estimate for K c {= w/ksT c ) and looking at the cu- 
mulant there (this is not shown here for brevity). In the 



thermodynamic limit we obtained U* = 0.57(1). This 
value is consistent with the cumulant crossings shown in 
Fig. HI 

In order to discard the possibility that the phase tran- 
sition is a first-order one, the energy cumulants [Eq. ((9])] 
have been measured. As it is well-known, the finite-size 
analysis of Ue is a simple and direct way to determine 
the order of a phase transitio n 6 ' 46 ' . Fig. O illustrates 
the energy cumulants plotted versus fc^T /w for differ- 
ent lattice sizes ranging between L = 24 and L = 60. 
The values of the parameters used in the MC runs were 
rn = 21, ni = 10 5 , n 2 = 10 5 , n M cs = 3 x 10 5 , T x = 0.25 
and T m = 0.15. As it is observed, Ue has the charac- 
teristic behavior of a continuous phase transition: the 
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minima in the cumulants tend to 2 /3 as the lattice size is 
increased. This indicates that the latent heat is zero in 
the thermodynamic limit, which reinforces the arguments 
given in the paragraphs above. 

Next, the critical exponents will be calculated. As 
stated in Rcfs. 48 49 50, the critical exponent v can be 
obtained by considering the scaling behavior of certain 
thermodynamic derivatives with respect to the inverse 
temperature K, for example, the derivative of the cumu- 
lant and the logarithmic derivatives of (8) and (S 2 ). In 
Fig. [BJa) we plot the maximum value of these derivatives 
as a function of system size on a log-log scaled. The 
results for l/v from these fits are given in the figure. 
Combining these three estimates we obtain v = 1.16(3) 
(see Table II). Once we know v, the critical exponent 
7 can be determined by scaling the maximum value of 
the susceptibilit y 48 ' 49 . Our data for xlmax are shown in 
Fig. E^b). The value obtained for 7 is indicated in the 
figure and listed in Table II. 

On the other hand, the standard way to extract the 
exponent ratio (3/v is to study the scaling behavior of 
(6) at the point of inflection, i. e., at the point where 
d(S)/dK is maximal. Since these points should scale as 

usual, (K$ - K c ) L 1 ^ = eL 1 ^ = const, we expect^ 9 - 

Wl lnf = i- /3/ ^(ei 1/,y )ocL-^, (14) 

where (5)| inf is the value of (8) at the point of inflec- 
tion. In addition, since the derivative with respect to K 
picks up a factor L x l v from the argument of the scaling 
function 8, 

= L (-/3/,+i/,)^/ (eL i/, ) K L a-P)/ r (15) 

dK max 

The scaling of (<5)| inf is shown in Fig. [6jb). The lin- 
ear fit through all data points gives /j(^'inf) = 0.24(6). 
In the case of d{8)/dK\ max [see Fig. Etb)], the value ob- 
tained from the fit is /3(<W/<"<W) = 0.20(8). Combining 
the two estimates, we obtain the final value f3 — 0.22(6), 
which is indicated in the figure and listed in Table II. 

The finite-size scaling theor y 13 i 36 i 37 i 48 allows for vari- 
ous efficient routes to estimate T c from MC data. One of 
this method, which was used in Fig.|H is from the temper- 
ature dependence of U (T) for different lattice sizes. An 
independent procedure to determine T c will be used in 
the following analysis. The method relies on the extrap- 
olation of the positions K C (L) of the maxima of various 
thermodynamic quantities, which scale with system size 

m- 13. 36. 37.48 

K C (L) = K c (oo) + const.L" 1/l/ . (16) 

Fig.[6jc) shows a plot of K C (V) vs. L~ x l v for the maxima 
of the slopes of (8), [d{8)/dK) m ^ U, (dU/dK) max ; 
ln<5>, (d\n(S)/dK) maK ; ln<5 2 ), (dln(S 2 ) /dK) max , as well 
as of the susceptibility, Xmax- The lines are fits of the 
data to Eq. (|16p with v = 1.16. From extrapolation one 



obtains i*Q (00) [or k B Tc (00) /w] for the different 
observables A. In this case, k B T^ U) '(oo) /w = 0.1824(12); 
k B T^ n{S) \oo)/w = 0.1825(14); fc B T c (ln<52>) (oc)/u> = 
0.1825(16); k B T^ 5) \oo)/w = 0.1820(20); and 

k B T^ x) '(oo) / 'w = 0.1823(5). Combining these esti- 
mates we find a final value k B T c /w — 0.1823(7), which 
coincides, within numerical errors, with the value 
calculated from the crossing of the cumulants. 

Strong corrections to scaling were observed to be 
present at small lattices and we excluded L < 60 data 
from the calculations. On the other hand, the L > 60 
data are not good enough to include corrections to scal- 
ing in the estimation of the critical exponents and critical 
temperature. Thus, the quoted errors in our results do 
not include systematic errors due to corrections to scaling 
that possibly could affect our data. 

The scaling behavior can be further tested by plotting 
(5)LP^ vs lelL 1 /", xL~ l/v vs eL 1 ^ and U vs |e|i 1/,y 
and looking for data collapsing. Using our best estimates 
k B T c /w = 0.182, v = 1.16, = 0.22 and 7 = 1.98, we 
obtain very satisfactory scaling as it is shown in Fig. [7] 
This study leads to independent controls and consistency 
checks of the values of all the critical exponents. The 
collapses in Fig. [7| were calculated by following a conven- 
tional FSS scheme. 

In the next, the analysis of Figs. [6] and [7] is repeated, 
this time applying the extended FSS scheme mentioned 
before. As it is shown in Fig. [51 the values obtained 
for the critical temperature k B T c /w — 0.1822(7) and the 
critical exponents v — 0.18(3), (3 = 0.23(6) and 7 = 
2.04(5) (see Table II) are in excellent agreement with 
those calculated using the conventional FSS. In addition, 
we apply the extended FSS scheme to collapse the data 
by plotting U , (8) and x m terms of the variables a = 
1-T C /T and Ly/T. The results are shown in Fig. U 
The behavior of the critical quantities for the extended 
FSS scheme is consistent with the behavior observed by 
following the conventional FSS scheme, which reinforces 
the robustness of the scaling analysis introduced here. 

The critical exponents obtained by conventional and 
extended FSS, along with the fixed point value of the 
cumulants, U* — 0.57(1), obtained in Fig. QJ suggest 
that the phase transition occurring for repulsive dimers 
on square lattices at 2/3 monolayer coverage belongs to 
a new universality class. We say "suggest" because the 
lattice sizes studied here do not allow us to exclude a 
more complex critical behavior, for example, the pres- 
ence of a tricritical poin t 34 ' 52 ! 53 . In order to elucidate 
this point, future studies on the whole phase diagram, 
including variables as coverage and an external field (cou- 
pled to the order parameter) breaking the orientational 
geometry of the phase, will be carried out. 

Finally, we will briefly refer to the specific heat expo- 
nent a. The "roughness" of the curves of C prevents a 
direct determination of a. However, the usual hyperscal- 
ing relations inequalities of Rushbrooke, a + 2/3 + 7 > 2, 
and Josephson, dv + a > 2 (being d the dimension of 
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TABLE II: Critical exponents obtained from the fits shown 
in Figs. E] and H 



exponent 


conventional 


extended 




FSS 


FSS 


V 


1.16(3) 


1.18(3) 


P 


0.22(6) 


0.23(6) 


7 


1.98(7) 


2.04(5) 


a 







the space), predict a negative specific heat exponent 
a w —0.3. This finding is consistent with the prelimi- 
nary results obtained in the present study. Namely, even 
though the fluctuations in the energy are of the order 
of the statistical errors in the simulation results (reason 
for which the data are not shown here), it is possible 
to observe that the maxima of the curves of the specific 
heat remain practically constant as L is increased. In 
other words, the specific heat seems not to diverge on 
approaching the transition. However, the present data 
do not allow us to be conclusive on this point and more 
work is needed to clearly elucidate the specific heats be- 
havior of the model. 



In summary, we have addressed the critical properties 
of repulsive dimers on two-dimensional square lattices at 
2/3 coverage. The results were obtained by using ex- 
change MC simulations and FSS theory. The choice of 
an adequate order parameter [as defined in Eq. ([2])] along 
with the exhaustive study of FSS presented here allow 
us 1) to confirm previous results in the literature ; 23 ! 24 
namely, the existence of a continuous phase transition 
at 2/3 coverage; 2) to calculate the critical temperature 
characterizing this transition; and 3) to obtain the com- 
plete set of static critical exponents for the reported tran- 
sition. Though it is not possible to exclude the existence 
of a more complex critical behavior, the results suggest 
that the phase transition belongs to a new universality 
class. 
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